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We develop a detailed analysis of synchrotron radiation in the effective Lorentz invariance violating 
(LIV) model of Myers-Pospelov, considering explicitly both the dynamics of the charge producing 
the radiation and the dynamics of the electromagnetic field itself. Within the radiation approxima- 
tion we compute exact expressions in the LIV parameters for the electric and magnetic fields, the 
angular distribution of the power spectrum, the total emitted power in the m-th harmonic and the 
polarization. We also perform expansions of the exact results in terms of the LIV parameters to 
identify the dominant effects, and study the main features of the high energy limit of the spectrum. 
A very interesting consequence is the appearance of rather unexpected and large amplifying fac- 
tors associated with the LIV effects, which go along with the usual contributions of the expansion 
parameter. This opens up the possibility of looking for astrophysical sources where these amplify- 
ing factors are important to further explore the constraints imposed upon the LIV parameters by 
synchrotron radiation measurements. We briefly sketch some phenomenological applications in the 
<t^ \ case of SNRs and GRBs. 

I/"") \ PACS numbers: 11.30.Cp,41.60.Ap, 03.50.Kk, 95.30.Gv 



in 
o 
o 

(N 



in 
o 
in 
o 

Oh 
i 

o 

> 

X 

5-H 



I. INTRODUCTION 
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The isotropy together with the Lorentz covariance of the laws of nature in locally inertial frames are properties 
which must be ultimately determined on experimental grounds. This is particularly true for boosts which, due to 
the non-compact character of the Lorentz group, can always be tested at higher energies. The experimental study 
of these questions began around 1960 with the works of Hughes et. al. and Drever p] and has continued since then 
involving a huge variety of very different testing methods with ever increasing experimental precision . It has been 
possible to correlate such highly disparate observations in terms of the Standard Model Extension (SME) proposed in 
\ Ref. H . This model is an effective theory that parameterizes all the possible Lorentz violating interactions consistent 
with the standard model, which arise in a spontaneously broken formulation of such symmetry violation motivated 
by string theory £| . The SME has been recently generalized to incorporate gravity effects [B| . 

A possible physical realization for such Lorentz invariance violations has recently been suggested within the context 
of dynamical modifications induced by quantum gravity in the propagation of particles in flat space, as for example in 
Ref. . Most theories of non-perturbative quantum gravity predict that our notion of space as a continuum manifold 
has to be abandoned at distances of the order of the Planck length £p w 1/Mp w 10 -33 cm, where space should be 
regarded as a bunch of quantum excitations. A first question related to this idea is whether or not such a drastic 
change in the description could produce observable residual modifications to our large scale description of particle 
dynamics, manifested as additional terms in the corresponding actions that are suppressed by powers of lp. Some 
approaches suggest that standard Lorentz covariance is consistent with the granular structure of space predicted by 
quantum gravity 0,01' expressed as the absence of ^-corrections to large scale dynamics, while others propose that 
such corrections should arise, but with different interpretations of their manifestations. The most direct interpretation 
of such corrections is through the existence of a preferred reference frame associated to the space granularity in which 
the existence of a minimum length provides a maximum physical upper bound for the momenta accessible in such 
a frame. Nevertheless, in Ref. |9j it was shown that this interpretation leads to severe fine-tuning problems when 
considering the radiative corrections induced by the standard interactions. Supersymmetry has been proposed as a 
way to improve such situation |Toj . An alternative proposal can be found in Ref. [Ll| . A second possibility arises when 
considering the dynamical modifications as signaling a spontaneous breaking of the Lorentz symmetry. In this case 
the corresponding effective theory also possesses a class of preferred frames that are those where the chosen vacuum 
expectation values of the respective fields are only a very small perturbation with respect to the Lorentz symmetric 
case. These are the so called concordant frames |12j and they are not directly related to the idea of a maximum 
momentum. The SME together with effective models like those of Refs. jToL Il3| belong to this category. A third 
possibility, generically included in the so called deformed special relativity approach, allows for dynamical corrections 
without the introduction of preferred reference frames by g enerali zing the concept of relative inertial frames on the 



basis of non- linear representations of the Poincare group |14lll5lll6lll7lll8j . Finally, alternative ways of incorporating 
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such dynamical modifications can be found in Ref. |19j . For recent reviews of the above topics see for example Ref. 
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[201 121|, |22J. 

One of the most studied consequences of the Lorentz violating corrections to the dynamics are related with multiple 
phenomena associated to the modified dispersion relations (MDR) for photons and fermions 

c, 2 (fc) = fc 2 ±^, E 2 (p)=p* + m 2 + r lR , L ^, (1) 

respectively. Here M is a scale that signals the onset of LIV and is usually associated with the Planck mass Mp. In 
this context the main emphasis has been set on the determination of observational bounds for the LIV parameters 
£ and r)R t L appearing in the MDR, thus opening the door to quantum gravity phenomenology j^. A partial list of 
references is given in p lIMIM I ^ . 

The possibility of deriving such modified dynamics from a more fundamental theory has been further explored in 
Loop Quantum Gravity, where constructions of increasing degree of sophistication have managed to produce modified 
effective actions for photons and fermions encoding the relations (JIJ . Basically, the starting point of these approaches 
has been the well defined Hamiltonian operators of the quantum theory, together with a heuristical characterization 
of the semiclassical ground state. The first derivation of such a type of dispersion relations in the context of Loop 
Quantum Gravity was developed in Ref. [2^| , and leads to a consistent extension of Maxwell electrodynamics with 
correction terms linear in lp. Subsequently, an alternative approach inspired by Thiemann's regularization p9| 
reproduced and extended such a result for the photon case |3(| and produced the corresponding modified Hamiltonian 
for the fermion case [IJ. For a review see for example Ref. [12 ■ Coherent states [33| have also been constructed, 
leading to similar types of dynamical modifications. All previous approximations have been kinematical, i.e. without 
incorporating the dynamical constraint of general relativity. In other words, a physically sound characterization of 
the semiclassical ground state corresponding to a given asymptotic classical metric that dictates the appropriate 
symmetries of the effective low energy theory is still lacking. This issue remains one of the most important open 
problems in loop quantum gravity. String th eory has also provided a possible connection between quantum gravity 
and LIV in the case of photons and fermions |34l l35| . 

Let us make some very general comments regarding some of the bounds already found in the literature for the 
LIV parameters appearing in MDR of the form given in Eq. Q for photons, electrons and protons, arising from 
dimension five corrections to standard electrodynamics. The main scenarios that we consider here are phenomena 
related to observations of electromagnetic radiation arising from astrophysical objects and from the arrival of ultra 
high energy cosmic rays (UHECR). 

The absence of birefringence effects from distant galaxies provides the bound |£| < 10 -4 j3||. Observations from 
the Crab nebula, a type of supernova remnant (SNR), give the following bounds on the parameters. On one hand the 
birefringence limit has been improved up to |£| J$ 10 -6 for the parameter region \rj e \ < 10 -4 with rj e being at least 
one of the t)l.r for electrons or positrons 22]. On the other hand, the intersection of the corresponding synchrotron, 
synchrotron-Cerenkov and the inverse Compton constraints produces the two-sided upper bound |?7l,jj| < 10~ 2 on 
one of the two parameters t}l,r 22] . Typical upper bounds for the corresponding energies in the Crab nebula are: 
10 6 GeV for the radiating electrons (positrons), 10™ 1 GeV for the synchrotron photons and 5 x 10 4 GeV for the 
inverse Compton photons. These give an idea of the energy scale in which the above bounds are obtained. 

The assumption of a dispersion relation of the type Q with a universal LIV parameter rj for all particles in the 
study of UHECR has produced the stringent bound \r}\ < 10~ 32 GeV -1 [37|. This reference also reports fits of the 
UHECR spectrum starting from 10 18 eU which include the AGASA data. 

The vacuum Cerenkov effect becomes significant for charged particles with high energy. It has been already 
discussed, both using a kinematical approach |26t l38l l39| , and on the basis of a dynamical analysis in the context of 
the Maxwell-Chern-Simons action [40(. An analogous analysis could be interesting in the Myers-Pospelov model, but it 
is out of the scope of the present work. Very stringent boundaries for the LIV parameters from the absence of vacuum 
Cerenkov radiation have been recently presented using protons in the ultra high energy region E prot0 n — 3 x 10 n GeU, 
based on reasonable assumptions about the nature of the high energy cosmic rays. Assuming that the values for the 
LIV parameter r\ are comparable for different standard model elementary species, or at least the same within species 
of a given spin, leads to consider unmodified dispersion relations for a structureless proton, i. e. r\ v = 0. With these 
assumptions the bound ||| ;$ 2a p = {^rn 2 proton / E^ roton ) ~ 10~ 34 GeV^ 1 is obtained in Ref. |33- Allowing for MDR for 

the proton, fj p ^ 0, the less stringent bounds of Ref. (2(| are recovered: fj p < a p /4, |£| < fj p — a p — a p ^l — Ar\ p ja 

Here f\ p = r\ p jMp and refers to the proton, while £ = £/il/p. These constraints apply to both photon helicities, which 
is the origin for the bound on the absolute value |£|. Considering the proton polarization we can write f] = X p k p , 
where A p = ±1 is the corresponding polarization. Thus, if we have k p > a p /4 (k p < — %>/4) the protons with positive 
(negative) hclicity emit vacuum Cerenkov radiation, while the protons with opposite helicity do not radiate. Hence, 
to establish boundaries from Cerenkov radiation it is necessary not only to know if there are protons at a given 
energy, but also what polarizations they have. At energies E ~ 10 11 GeV the UHECR are detected by the cascade 
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they produce in the atmosphere, where the information about primary proton polarization is lost. To determine their 
polarization they must be observed directly as primary cosmic rays. The highest energy for which this could be done 
is E ~ 10 5 GeV, with instruments on board of satellites |4lj |. In this case a p ~ 10 -15 GeV -1 and thus we would 
get |£| ;$ 10 , which is 10 10 times weaker than the boundary obtained from vacuum birefringence. It is interesting 
to remark that the highest electron energies measured in primary cosmic rays are estimated to be of the order of 
10 3 GeV with a flux 0.05 times that of the protons j^J. This gives a e = 10" 15 GeV -1 . In the Crab nebula we have 
a e — 10 -24 GeV -1 which produces a much better bound than electron vacuum Cerenkov radiation in cosmic rays. 

For a detailed and updated discussion of the LIV parameter constraints see for example Refs. |20Ll2lLl22| . 

Now let us go back to the observation of 100 MeV sinchrotron radiation from the Crab nebula. This fact alone 
leads to the bound ±t}l,r > —7 x 10~ 8 for at least one of the four possibilities [25j. By itself, this does not impose 
any constraints on the parameters, but it plays a role when combining it with the vacuum Cerenkov effect. Such a 
bound is based on a set of very reasonable assumptions on how some of the standard results of synchrotron radiation 
extend to the Lorentz non-invariant situation. This certainly implies some dynamical assumptions, besides the purely 
kinematical ones embodied in p. In this paper we examine these assumptions in the light of a particular model, 
which we choose to be the classical version of the Myers-Pospelov (MP) effective theory that parameterizes LIV 
using dimension five operators 0. Furthermore, a complete calculation of synchrotron radiation in the context 
of this model is presented. This constitutes an interesting problem on its own whose resolution will subsequently 
allow the use of additional observational information to put bounds upon the correction parameters. For example 
we have in mind the polarization measurements from cosmological sources. The case of gamma ray bursts (GRB) 
has recently become increasingly relevant although it is still at a controversial sta ge Epf - Our calculations rest 
heavily upon the work by Schwinger et al. on synchrotron radiation, reported in Refs. |44l l45l |46| . In this work we 
restrict ourselves to classical theories, and therefore we do not address the fine-tuning problems arising from quantum 
corrections associated with the coexistence of space granularity and preferred reference frames 0- A partial list of 
previous studies of electrodynamics incorporating LIV via dimension three and four operators is given in Refs. [40ll47|. 
A summary of the present work has already been reported in Ref . |48j . 

The paper is organized as follows. In Section II we obtain the equations of motion for the charge and electromagnetic 
sectors of the MP model. The polarized electromagnetic fields for arbitrary sources in the radiation approximation 
are subsequently calculated in Section III using the standard Green function approach. In Section IV we obtain 
the general expression for the angular distribution of the power spectrum. This is subsequently applied in the next 
section to the case of the synchrotron radiation produced by a charged particle moving in a plane perpendicular 
to a constant magnetic field. The corresponding angular distribution of the radiated power in each harmonic is 
obtained, together with the corresponding total power. Also the associated Stokes parameters are calculated in this 
section. The expansions of such exact results to leading order in the LIV electromagnetic parameter are contained 
in Section VI, where we highlight the dominating effects of the Lorentz violation. A rough estimation of the relative 
contributions arising from dimension six operators is included in seccion VII. Motivated by the estimated parameters 
of the astrophysical objects under consideration we discuss the large harmonic expansion of the previously found 
results in Section VIII. We analyze some phenomenological consequences of our results in Section IX. Finally, Section 
X contains a summary of the main results obtained in the paper together with a discussion of accessible regions 
for the LIV parameters in the cases of Crab Nebula, Mk 501, and the 7-ray burst GRB021206. There are also two 
Appendices containing frequently used material. 

II. THE MYERS-POSPELOV ELECTRODYNAMICS 

The Myers-Pospelov approach is based on an effective field theory that describes Lorentz violations generated 
by dimension five operators, parameterized by the velocity of a preferred reference frame which is not taken as a 
dynamical field. These operators are assumed to be suppressed by a factor Mp 1 , and hence can be considered as small 
perturbations at the classical level. In the following we analyze the electromagnetic radiation of a classical charged 
particle in this framework. As a first step in this approach we characterize the particle and the electromagnetic field 
dynamics, introducing the charge-electromagnetic field interaction via the minimal coupling, in such a way that the 
usual gauge symmetry is maintained in this theory. 
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A. Charged particle dynamics 



The dynamics of a classical charged particle can be obtained from the action for a scalar charged field. In this case 
the Myers-Pospelov action is 



Smp —Id 



d^d^tp - n 2 (p*(p + ifjip* {V ■ df ip 



(2) 



with the notation V ■ d — d^. In momentum space, where we write f(x) = ipo ex.pi(p°t — p • x), and in the 
reference frame where V a — (1,0), the dispersion relation becomes 

(p°) 2 +v(p°) 3 = P 2 +» 2 - (3) 
To compare with Jacobson et al. results l2(jL |2?J , the parameter rj can be written 



'/ 



(4) 



where Mp is the Planck mass and rj is a dimensionless constant. The equation is an exact relation in fj. From 
here we obtain the Hamiltonian for a massive particle to second order in fj 



P ° = ff=(p 2 + A1 2 ) 1/2 -ir;(p 2 + / , 2 ) 



.3/2 



o(r) 



(•5) 



In the following we will consider the interaction of a particle of mass fi and charge q with a static magnetic field 
(4> = 0, A(r)). The standard minimal coupling produces the Hamiltonian 



H 



1/2 



1 



t>~2 
o 7 ? 



c 



3/2 



(6) 



Here c = 3 x 10 10 cms 1 denotes the uncorrected velocity of light in vacuum. In the sequel we set c = 1. Observe 
that the dispersion relation © provides the exact inversion 



(p-qAf = (l + fiE) E 2 -^ 2 , 
with E being the energy of the particle, which together with the Hamilton equations yield 



dH _ 1 

Wi ~ E 
dH _ . dAj 

dxi dxi ' 



3 9 

*< = v i = ^Er = 9 ( 1 - + I 7 ' 2 ^ 2 ) ^ ~ ^ ! 



Pi 



up to second order in fj. The acceleration results 
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r = 



E 



l-^E+^-fE 2 ) (v B ) 



(7) 

(8) 
(9) 

(10) 



As in the usual case, this means that the magnitude |v| is constant and that the projection of the orbit in a plane 
orthogonal to B = V x A is a circular orbit with a Larmor frequency 



UJ = 



\q\B 

E 



1 -lfjE+-rj 2 E 2 
2 ' 4 ' 



(11) 



In general the motion is a helix with pitch angle (the angle between the velocity and the magnetic field) a. Let us 
define 



[3±_ = /3sina, /3\\ = /3cosa, 



(12) 



where we emphasize that we are using the standard definition (3 — |v|/c. The solution to the equations of motion can 
be written as 



r(t)= ( — coscjoi, -^sinwot, (3\\t) , 
\uj ujo J 



(13) 
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and hence 



v(t)=(— (3± sinuJot, f3±cosuj Q t, /3m). 



The orbit equation l|13|) identifies 



iv 



as the radius of the projection of the helix in the plane perpendicular to B. From Eqs. J7|) and © we obtain 



P E 2 



1 



f)E 3 15 rfE 4 



where the considered range of energies is such that 

/' 



E 



« 1, f]E « 1. 



(14) 



(15) 



(16) 



(17) 



According to the relation among the quantities appearing in Eq. (|17|) . we can consider different energy ranges to 
write approximate expressions for the Lorentz factor 



7 = (l ~ ft 2 ) 



2\-V2 



(18) 



In Eq. 1)16(1 . the ratio between the second order term in fj and the first order one is proportional to fj E ^ 1. Therefore, 
assuming fj > 10 -19 GeV -1 , the last term is negligible compared with the second one at least for E < 10 19 GeV. The 
last term is of order one when ry 2 i? 4 //! 2 ~ 1, i.e. when E ~ (/i/77) 1 ^ 2 . For electrons this means E ~ 7 x 10 7 GeV ~ 
7 x 10 4 TeV, and an even larger energy for more massive particles. Therefore, in the range of energies expected in 
astrophysical objects and when E < 7 x 10 7 GeV, the 7 factor can be written as 



7 



E 



1 + 2^) 



-1/2 



Hence we have the two possible approximate expressions 

fjE 3 



7 



7 



A* 



— fi— ^ 

2k£ V ^e 3 J 



< 1, 



> 1. 



(19) 

(20) 
(21) 



The limiting condition between these two energy ranges is E ~ 10 4 GeV for electrons and £ ~ 10 6 GeV for protons. 
In the case of electrons in a SNR of the type of Crab Nebula, the maximum energy of the electrons is of the order of 
l0 3 TeV, and therefore the approximation l|20|l still holds. 

According to the preceding analysis, the current for a charged particle moving in the general helical motion is 



j(t,r) =qS 3 (r-r(t))v(t), 



(22) 



where r(t) and v(t) are given in Eqs. (|13fl and l|14fl respectively, with /3 determined by the Lorentz factor (|20|l . In 
the following we will consider only circular motion, i.e. a — n/2. 



B. The electromagnetic field 

The Myers-Pospelov action for the electromagnetic field is 



Smp = J d A x —F^F"" - 4Kj»A„ + £(y a F ai ) (V ■ d)(V j3 F l3S ) 



where £ = £/Mp, with £ being a dimensionless parameter. We are using the conventions 
F^ = \e^F pa , e 0123 = +l, e 123 = +l, f) = (+,-,-, -), 



(23) 



(24) 
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As usual we have = d^Ay — d v A^ and j M = (p,j), with the electric and magnetic fields identified according to 

Foi = Ei, Fij — —ajkBk, Bk = —-^tkijFij, (25) 



Foi — Bi, F^ = eijkEk, Ef. — -CkijFij. (26) 



These lead to the standard homogeneous Maxwell equations 



V-B = 0, VxE+ — =0. (27) 

The equations of motion obtained from the action 1|23|1 are 

d^ v + 1{V ■ d) (e^Vp (V a d p F aS ) - (V a d a ) {VpF?")) = 4nrj v . (28) 

For simplicity, we keep working in the rest frame V a = (1, 0). If required we can boost to any other reference frame 

with a given V a by means of an observer (passive) Lorentz transformation. Rewriting Eqs. I|28|) in vector form we 
get 

V ■ E = 47rp, (29) 

? + VxB + 4f-VxE + ^W (30) 



dt "dt\ dt , 

Thus, in terms of the standard potential fields — (<fr, A) we are left with 

„o, d 



, , gt V-A = -4wp, (31) 
d ( dA \ d 2 A 

aU w + ^) +Vx(VxA) + 2eVx ^ r = 4?rj ' (32) 



which can be conveniently expressed in the radiation (Coulomb) gauge, V • A = 0, as 



1 

V 

d 2 A 2 . . „ r „ 8 2 A . ( . „ 1 



-4tt^2P, (33) 



dt 2 s dt 2 V V 2 

where the electric and magnetic fields reduce to 

dA 



V 2 A + 2£Vx— r =4^ j-V-^V-j =4^j T , (34) 



E = - — , B = VxA, (35) 



in the corresponding large distance approximation. 

The energy momentum tensor T^ v for this modified electrodynamics is given by 



S = -ExB-fEx- (37) 
4?r 4tt dt ' v ; 

which are exact expressions in £. These components satisfy the usual conservation equation 

g T 00 



dt 



V • S = 0, (38) 



outside the sources. 

To solve the equation of motion for A, Eq. I|34|) . we can go to the momentum space with the convention 

F(t,r)= ( -^Le-^+^F^k). (39) 
J (2,1:) 
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The different types of Fourier transforms are specified by the corresponding arguments. For example, if F(t, r) denotes 
the function in space-time, F(u>, r) denotes the Fourier transformed function to frequency space, while F(u>, k) denotes 
the Fourier transformed function to frequency and momentum spaces. In this way Eq. ij^Hl reduces to 

(-uj 2 + k 2 - 2ilu? kxj A(w, k) = 47rj T (w, k). (40) 

This equation can be diagonalized using the basis of circular polarization A = ±1, defined in Appendix A, giving 

(-u? + k 2 T 2£w 2 fc) A ± (w, k) = 47rj±(w, k). (41) 



III. GREEN FUNCTIONS AND FIELDS 



The simplest way to proceed is by introducing the retarded Green functions with definite polarization 

1 



G£ et (w,k) = 



k 2 - \2£u 2 k - lo 2 



(42) 



together with the total retarded Green function [G re t(to, k)]ik = J2\ ^ik Wet( w i k) an d to calculate [G re t{uj, r — r")]ifc. 
Here Ph. is the helicity projector defined in Eq. (|A4|) . The angular integration can be immediately performed and the 
remaining integration over dk gets contributions only from the upper half plane in the complex variable k. Finally 
one can identify the polarization components of the total Green function as 



J (27T) 



-k) 



1 



n(Xz) 



in(Xz)ujB. 



AirRy/TTz 2 

where R = |r — r'|. Here we introduce the polarization-dependent refraction index n(Xz) 



i(Xz) — yl + z 2 + Xz, z = £u 



(43) 



(44) 



In this way, the fields A A in Eq. (|41|l have well defined phase velocities v\ = l/n(Xz) and this situation can be 
described as the propagation of photons in a dispersive birefringent medium. 

The Green functions 1)43(1 determine the corresponding potentials with the standard replacements 1/R ~ l/|r| = 1/r 
in the denominator and R ~ r — fi • r' in the phase, with n = r/r being the direction of observation. Note that in the 
Green function phase 



n(Xz)u> |r 



cur 



n . r - n • r 

1 h X£u - X£uj 



(45) 



we are taking the radiation approximation, which means that the subdominant terms of order (r'/r) 2 or higher are 
neglected. Consistency demands that terms proportional to the LIV parameter £ are larger than the neglected one in 
order to properly include them in this phase. Our general results are presented in this full far-field approximation. 
Using iffify we obtain 

A A (cj,r) = 4vr / riV G^ t (uj, r - r') r') e m(Xz ^ r f dV e -<[n(A*)aA]V r ') ; ( 46 ) 

J ' R\/l + z 2 J 



and thus we finally get 



A>,r) = l -^M=e m ^]\^ x ) 

r VI + z 2 



(47) 



in the radiation approximation. The fields A + (w,r) and A (u>,r) correspond to right and left circular polarization 
respectively. Let us emphasize that the momenta 



\i\ = n{Xz) lu h 



(48) 



in Eq. (|47|l are fixed in terms of the frequency and the direction of observation. As usual, the integration over 
has been conveniently written as the Fourier transform in momentum space jj,{u!,kx) of the function jj,(w,r') in 
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coordinate space. In the RHS of Eq. I|47|l we have used the relation (|A1|) for the transverse current. The full vector 
potential is given by the superposition 



A(w,r) 



1 1 



n(Xz)e m ^ r i x (cj,k x ) 



A=±l 



Hence the electric and magnetic fields are 



B(w,r) = V x A(w,r) 

_ 1 UJ 

E(w, r) = iwA(u, r) 

1 ILU 



2 (z)e™ (z)wr j+(w,k + ) - n 2 (-z)e™^ z ^ r j-(w,k_) 



l (z)e m ^j+(^k + ) + n(-z) e m (- 2 ^j-^,k_) 



r VTTz 2 

Note that, contrary to the standard case, where n x E oc B, here we have 

1 



n x E(w, r) 



y/T+z* 



[B(w,r) + izE(w,r)] . 



(49) 



(50) 
(51) 

(52) 



IV. THE ANGULAR DISTRIBUTION OF THE POWER SPECTRUM 

This is defined as , \J , where P(T) is the radiated power at time T into the solid angle dSl in a given radiation 
problem. We can compute the total energy emitted in terms of the Poynting vector l|37|) 

/oo 
dtn-S(t,r) r 2 dn. (53) 
-CO 

This last expression can be rewritten introducing the Fourier transform of the Poynting vector, 

d 2 E 



E = 



f f d E f°° dui 

J duj J dn dnd^ = J 2^ [n ' s(cj ' r) + n ' s(_w ' r)] r2dQl 



and allows us to obtain the angular distribution of the energy spectrum 



d 2 E 



[n.S(«,r)+n-S(-«,r) 



dVLduj 2n 

from where the angular distribution of the power spectrum can be identified as 

d 2 E 



+ °° dT d 2 P(T) 



dfldui J -oo dwdVt 
The Poynting vector is given by Eq. H37[). from which results 

S(u, r) = — (B(-u, r) x B(u, r)+ioj£,E(-uj, r) x E(w, r) ' 

47T V J 

This expression, exact in £, can also be written in terms of the potentials, in which case we have 

S(w,r)= — — \uj 2 A{— uj, r) x n(£oj)A + (uj, r) — n(— |w)A_(w,r) 

47T L 

-|w 3 A(-w,r)xA(a;,r)| . 
Finally, using the relations (|A3(I of Appendix A we get 



(54) 



(55) 



(56) 



(57) 



(58) 



n- S(o;,r) = — \fl + z 2 [A_(-cj,r) ■ A+(w, r) + A+(-w, r) ■ A_(w,r) 

47T 



(59) 
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Introducing Eq. (|59() in Eq. I|55|) we finally have 
d 2 E r 2 u 2 



dfldiu Ait 2 



y/l + z 2 [A_(-w, r) ■ A+(w, r) + A_(w, r) • A+(-w, r)] . (60) 



Now we need to express the products A T (— r) • A±(w, r) in terms of the current j(w, k) via the relation l|47|) . To 
this end it is convenient to introduce the projectors (|A4|) that satisfy the relations (|A5|) . I|A6() and (|A7|) . Using these 
results together with the general relation 

&(-w,-k)=#(w,k), (61) 

we obtain 



A r) • A+ (±w, r) = 1 (w, um(±z)fi)) P± j,. U. ( ± : in i . ( (i2 ) 

which leads to 



= \h\Tz~ 2 [ n2[z) Jt (tJ ' k+) Pl+k M^+) + n2 (- z ) fk(^-) Pu 3i (w,k_) ] . (63) 

In order to identify the angular distribution of the power spectrum we need to go back to Eq. (|56|l . Each contribution 
in Eq. I|63l) is of the type 

C H = j* k (u, k) X kr j r («, k) = / dt dt'er 1 ^-' )j* k (t, k) X kr Jr (t\ k) . (64) 

J — oc 

Introducing the new time variables r — t — t 1 and T = (t + t f ) /2 we get 

/+oo /»oo 
<£T / dre^ T ^(T + T/2,k) X fcr jV(T-r/2,k). (65) 
-OO J —OQ 



Inserting this last relation in Eq. (|63|l and comparing with Eq. (|56() we obtain the final expression for the angular 
distribution of the radiated power spectrum 

J2p/-Ti^ i 2 poo 

+n 2 (- 2 )inr + r/2,k_) P t - r i r (r-T/2,k_) ] , (66) 
as the sum of the contributions of both circular polarizations. 



V. SYNCHROTRON RADIATION 



Let us rewrite Eq. (|66|) in the form 

d 2 P(T) _ ^ lu 2 n 2 (\z) 



.2 2M "i poa 

d^dn ^7^L dTe ~^ jUT + T/2M) P *ririT-r/2M), (67) 

where the current associated to the motion of a particle in an external magnetic field was already obtained in Eq. 
(J22J- Its spatial Fourier transform is 

j(t,k) = J d 3 y(t,r)e-* y = qv^e-^M , (68) 

and thus we get 

X=±l v ' 
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where 



/oo 
dre-^jt (T + r/2, k A ) Pi jr (T - r/2, k A ) . (70) 
-OO 

This integral can be expressed as 

J(T,CJ,X,n) = q 2 dre -«,r e «,„(A,)n.(r(T+r/2)-r(T-r/2)) 5) 



with 

= ^(T + r/2) P£ r v r (T-T/2) 
= - [v(T + r/2) • v(T - r/2) - n • v(T + r/2) n • v(T - r/2) - iAn • v(T + r/2) x v(T - r/2)] . (72) 

Here it is useful to rewrite the terms containing n • j « n ■ v in terms of the charge density. To this end we use 
current conservation in momentum space 

wp(w,k)-k.j(o; 1 k) = Q, (73) 

obtaining 

n.j(w,k A ) = —f — p{ui,k x ), (74) 
n(Xz) 



which leads to 



a.j(t,k A ) = -^- T p(t,k A ). (75) 

n(Xz) 



In this way we have 

1 

n 2 {Xz) 

In other words, in the expression of 3, Eq. l|72|l . we can make the replacement n • v(t) — > l/n(Xz), leading to 



n • j* (T + r/2, k A )] [n • j (T — r/2, k A )] = -3-77^* I T + T A k A] P [T - r/2, k A ] . (76) 



3 = - [v (T + r/2) • v (T — r/2) - n~ 2 (Az) — iAn • v (T + r/2) x v (T - r/2)] . (77) 

In the following we will consider the simplest case, when the orbit of the charged particle is in a plane perpendicular 
to the magnetic field, i.e. the case where /3y = and 0± — (3. To continue with the calculation we introduce the 
direction of observation 

n= (sin^ 0, cos 0), (78) 

in the same coordinate system where r(t) and v(t) were defined in Eqs. (|13fl and i|14|) . Recalling the corresponding 
expressions we calculate 

v(T + r/2) • v(T - r/2) = (3 2 cos u t, (79) 

n- v(T + r/2) x v(T - r/2) = ~(3 2 cos sin w r, (80) 

and substituting in Eq. I|77|> we are left with 

= - [/3 2 coscjot — n~ 2 {Xz) + iX(3 2 cos sincj ' r ] ■ (81) 

Now we can compute the cycle average over the macroscopic time T, of the angular distribution of the radiated 
power spectrum 



d 2 P{T) \ _ Lo 2 n 2 {Xz) 

dtudfl / ~ 4tt 2 VTT^ 



' A =±i 
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The dependence on T is in the exponential factor of Eq. (|71|l , where we have 

n- [t(T + t/2) -t(T-t/2)} = -2#sm0sin (w T) sin (co t/2) . (83) 
Thus the required averaged value is 

< e «,n(A,)fl.(r(T + r/ 2 )-r(T-r/2)) >= j^YsinO), (84) 

with 

Y = u n(Xz) R sin(w r/2). (85) 

From here we get 

2 /*oo 

<J{T,u,X,n)> = \ I dr e-^ T J (2Ysm6) 

^ J -oo 

x \p 2 cos uj t - n~ 2 (Xz) + iXj3 2 cos 8 sin (u t)] . (86) 
Now we apply the addition theorem for Bessel functions |4flj 

oo 

Jo (k VP 2 +° 2 - 2pacos(q> - <//)) = Jra{kp)J m (kcj)e ml{ ^'\ (87) 

m=— oo 

with 

p = a = R 7 k = ujn(Xz) sin 8, (0 —(/>')= lj t, (88) 

which yields 

oo 

Jo(2rsin6»)= ^ J^[wn(A2)i?sin6i]e irnwoT . (89) 



The next step is to rewrite the trigonometric functions depending upon uiqt in Eq. I|81|) as exponentials and to perform 
the r integration. In such a way we finally arrive at the expression for the observed time average of the power angular 
distribution, given by 

d 2 P(T) \ f A , /P m , A ^ Qn . 

' A=±lm=0 

with 

d^A = g^g rv^0_ | ^ [(1 + Xcqs6) j2 m _ i{Wmx) + (1 _ Acos0) J^ +1 ( WmA )] 

-2[n(Az m )]- 2 J^ (W mA ) } , (91) 
where o; TO = mujQ, z m = £,u) m and 

W m \ = uj m n(Xz m )R sin 8. (92) 

Using the relation loqR = (3 we can rewrite VF m A as 

V^ m A = mn(Xz m )/3 sin9. (93) 

We notice that in the limit £ = the unpolarized angular distribution of the power obtained from Eq. (|91|l coincides 
with the standard result given in Eq. (38.37) of Ref. [iEf. Let us observe that to first order in the LIV parameters 
we can take 

^o = e^ (94) 
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according to Eqs. (jnj and lO, since £fj is of second order in the correction parameters. 

Using additional properties of the Bessel functions and following Ref. 0] we can recast Eq. I|91|) as 



dP rn \ _ U m q n 2 (Xz m ) 



-, [ m [ {[PJ' m (W m x)f + [n-\Xz m )cot6 J m {W mX )} 2 \ 



n(Xz m ) iff> dJl(W mX ) 
4tt VI + z m 2 dW mX 
It is interesting to remark that the above expression can be written as a perfect square 



dP mX u 2 m q 2 1 



[X/3n(Xz m )JU(W mX ) + cot9 J m {W mX )} 2 . (96) 



dQ 4tt ^/TTzl 



A. Total spectral distribution 

To calculate the integrated power spectrum it is convenient to start from Eqs. H82[l and l|86(l . Integrating over the 
solid angle £1 we have 

where the required integral is 

[ sin 6d6 < J(T, u, A, n) > = %- [ dr e~ WT 

J * J-oo 

x { [P 2 cos uj n T - rT 2 (Xz)] h + X0 2 (i sin uj t) I 2 } . (98) 

Here 

r 71 sin 2Y 

h= dd sine J o (2Fsm0) = — — , (99) 

.in y 



h = [ dO sin 9 cos J (2Y sin 6) = 0, 
Jo 



(100) 



where Y is given in Eq. (|85|) . The last integral, related to parity violation, is zero by symmetry Substituting in Eq. 
(|97jl we get 

where 

= sin[2^n(A,) sin(0/2)] ^ _ n _ 
sm(0/2) 

Following the method of Ref. |45| we introduce the Fourier decomposition of the periodic function f x (4>) 

h{4>) = J e^f Xm , f Xm = r g e -^/ A (0), (103) 

m— — oo ~ 71 

and performing the r integration we have 

/ ' dP(T) \ q 2 uj V*je/ , mn(Az) 

* ' A=±lm=0 v 

It only remains to compute f Xm 

t r d<t> -»mtf, sin[2cjTi(Az)i? sin (0/2)] f 2 , _a n M /, ncn 
hm = J_^ e sinW^ cos - n (Az)j . (105) 
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The only contribution from the exponential comes from the symmetric part under 
terms of the (f>/2 we are left with 



r* dd> , sin \2um(\z)R sin (6/2)] r ^ 9 9/> . „„ 9 . 9 , , 

/a™ = / — cos m0 ^ \ / VW yj /? 2 - n- 2 (Xz - 2/3 2 sin 2 0/2 . 

Jo sin (0/2) 

Using the following representations of the Bessel functions [ifij 

J2m(%) = / — cos 7710 cos (a; sin (6/2)) , 

JO T 

J' 2m (x) = — / — sin (0/2) cos 7710 sin (x sin (0/2)) , 



After rewriting cos0 in 
(106) 



dyJ2m(y) 



sin (a; sin (0/2)) 



7r sin (0/2) ' 

with x = 2Rum(Xz) = 2mRu>an(Xz), and recalling that z = £lu = ^mujQ and Rujq = (3, we arrive at 



f Xm = 2l3 2 J' 2m (2m[3 n(Xz)) - [ n - 2 (\z) ~ (3 2 



2mn{\z)P 



d,X J 2 m(x)- 



Thus the total averaged and integrated power radiated in the m harmonic is 



q 2 mu>o 



2RV1 + z 2 



jn(z) 
+n(—z) 



2p 2 J' 2m (2mpn(z)) - \ n - 2 (z) - /3 2 



2mn(z){3 



dx J2m(x) 



2p 2 4 m (2mPn(-z)) - \n- 2 (-z) - (i 2 



2mn( — z)j3 



dx Jl m (x) 



(107) 
(108) 
(109) 

(110) 



(111) 



This result is exact in z — £oj. Let us observe that each separate contribution in Eq. (|lllfl has the form given in Eq. 
(2.22b) of Ref. ^(|, after the appropriate change e 2 /47r — > q 2 is made, multiplied by the additional factor 



1 n(Xz) 



n(Xz) 



n 2 (Xz) 



2 V TTi 2 n(z)+n(-z) l + n 2 (Xz)' 



(112) 



respectively. The calculation performed in this reference is for a constant index of refraction but, as mentioned in 
footnote 3, the corresponding results hold for an u> and H dependence of n and fi. At this stage, the additional factors 
(|112fl together with the specific form of the refraction indices are the remainders that we are considering a different 
electrodynamics. 



B. Polarization of the synchrotron radiation 



In general, the electric field (|51|) . which is naturally written in terms of the circular polarization basis e± introduced 
in Appendix B, will exhibit elliptic polarization. To describe the polarization we will use the Stokes parameters 
according to the definitions of Ref. [51| 



/=|(e;.E(c,r))| 2 + |(e*_.E(c,r))| 2 
UH(e;.E(c,r))| 2 HK.E(a,,r))| 2 
Q = 2Re \(e% • E(w,r))* (e*_ • E(w,r)) 



U = 2/777 



«-E(w m ,r))*(ei-E(u; m ,r)) 



(113) 
(114) 

(115) 
(116) 



It is convenient to define the reduced Stokes parameters by dividing each one by the intensity /, in such a way that 
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satisfying the relation 

l = v 2 +q 2 +u 2 . (118) 

As usual, purely circular polarization is described by v = 1, q — u — 0, while purely linear polarization corresponds 
to v = 0. 

A first step in the description of the polarization is the explicit calculation of the electric field, which is given in 
terms of the currents k+) and k_). We need to calculate the time Fourier transform of the current 

for each polarization A and subsequently project into the circular basis. We start from 

j (t, k x )=qv{t)e- ikx - rW =qf3 (- sinwo*, cosw t, 0) e^"^, (119) 

where we can write 

e — ikx-r(t) _ e -i^n(\z)0 sin 9 cos u t (120) 

The Fourier transform is 

/oo 
dt j(t,k x ) 
-00 

= q/3 dt e luJt {-smoj t,cosuj t,0)e~ lWxCOS ' AJot , (121) 



with W\ = (u> n(\z)0 sin 0) hop. Using the generating function for the Bessel functions of integer order (see for 
example Eq. (38.27) in Ref. gg) we obtain 



g—iW\ cosujo* 



{-i) m e- im " ot Jm{W x ), (122) 

m— — oo 

and in this way we have 

j(w,k A )=g/3 Y, H) m Jm(Wx) e^- m ^(-smLu Q t, coscu Q t, 0) dt. (123) 

ro=— oo J -co 

Now we can perform the integration with respect to t. The final result can be written in terms of each mode 
contributing to the current 

oo 

j( w ,k A ) = Y 5 ( w_ma; o)jm(w m) k A ), (124) 

where 



m— — oo 



j m (w m ,k A ) = 2ir(~i) m ql3 



in -J m (W m \), i-jj- — J m (W m \),Q 



(125) 



with W m \ given by Eq. (JS^J ■ 

The projection j^j(w m ,k A ) of the currents into the circular basis is given by Eqs. IjAlOl! of the Appendix. Using 
these relations together with the definitions (|51|) we arrive at the final expressions for the circularly polarized electric 
fields 

E+(c m ,k + ) = -^ { -i)^A^A^ {Zm ). m r (B + + A + Q) e+j (126) 
E-( Wm> k_) = -V2nHr^^t^e"^^ (B m - A m co S 0) e_, (127) 

~ + Z m 

where the notation is 

j m (u; m ,k A ) =2n{-i) m q [A x m , iBi, 0] , (128) 
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with 



A ™~PWr — J m{VV m x), B m -P — -7TT- ■ 

W\ m aW m X 



The combination A cos 9 reduces to 



At cos = D 



x _ J m {W mX ) 
tan 6* 



(129) 



(130) 



In the standard situation, where £ = 0, the above expressions produce the result 



sm f 



Jm(W m \), W 



dJ m (W m \) 



dW, 



riiX 



, - COS 9J m (W m A) 



(131) 



which is clearly orthogonal to n and parallel to the corresponding electric field given in Eq. (6.65) of Ref. [53 for the 

case (3\\ = 0. 

Next we calculate the Stokes parameters, Eqs. 111161 . starting from Eqs. I|127fl . The required projections are 



« ■ E m ( Wm> r)) = -V2^r q -^4e™^ {B + + D + } 



r X /T+ z- t 

(e*_ • E m ( Wm ,r)) = -V2n(-irS. ^(- Z ^ e in { - Zm ).r (B - _ D - 



(132) 
(133) 



and we also need 



(e; • E m K„,r))* (el ■ E m (c m ,r)) = 2 "M^ e^M^K 



r J 1-1- 2^ 



(134) 



In this way we arrive at the reduced combinations 



l-Rl 



<ln 



2R r , 



)| [n(-z m ) - n(z m )] to m r\ 



2R n 



which completely characterize the polarization of the wave in terms of the polarization index 

n{-z m ) (B~ - D m ) 



in | [n(-z m ) - n(z m )] w m r| , (135) 



(136) 



From the above relations the verification of the condition 1118|) is direct. The purely circular polarization corresponds 
to R m = 0, while the purely linear polarization is given by R 2 n = 1. Substituting the corresponding expressions for 
B\\\ and D± we obtain the general expression for the polarization index 



Rrt 



n(-z m ) (pam6J^(W m -)-wa6J m (W m -)) 



n(z m ) {(3 sin 6 J^(W m +) + cos 6J m {W m+ )) 
In the limit £ = we have B^ = B m and Df^ = D m and Eq. (|136fl reduces to 

B m — D 

n 



R„ 



1 _ 12m. 

B m 



Bra ~\~ Dm 



1 



The additional definition 



leads to 



A, 

Br, 



J m (m/3 sin 6) 
(3 tan # J' m (m/3 sin 6*) 



2r„ 



1 + r 2 

1 ' n 



1 - r 2 

_ ' r) 

1 + r 2 , 



0, 



(137) 



(138) 



(139) 



(140) 



which coincides, up to the sign of r m , with the parameters introduced in Ref. [5 
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VI. DOMINATING EFFECTS OF THE LORENTZ VIOLATION IN PHOTON DYNAMICS 

Several characteristics of the synchrotron radiation are affected by the Lorentz violating parameter £. From the phe- 
nomenological point of view the more relevant are the power angular distribution, the total power and the polarization 
of the radiation. In what follows we discuss these effects. 



A. Angular distribution of the power radiated in the unpolarized m -harmonic 

In this subsection we expand the distribution l|95|) to first order in £. Notice that its general structure is given by 

dP m ,\ 



F m (AO+AG m (AO, 



and thus we get 



dP„ 



F m (0) + £G' m (0) + A |>4(0)+G m (0) 



(141) 



(142) 



where the derivatives are with respect to the argument A£ of the functions. When we sum over polarizations it is 
evident that the term linear in A vanishes, but nevertheless there remains a contribution proportional to £. Here we 
identify the unpolarized average angular distribution of the power radiated into the m th harmonic 



dP r , 



dP, 



dP„ 



= 2 



F m (0) + £G' m (0) 



dn dn dn 

From Eq. (|95(l . and eliminating the second derivative of the Bessel function via its differential equation we get 



F m (0) 



47T 



where 



{[pj' m {w m f + [C0t9 J m (w m )] 2 } , 

= (q^rmwm e r j L{Wm)] 2 + (1 _ p sin 2 , 

Z7T | 

W m = W m A|f_o = mfisme. 



J m (W m ) 



(143) 

(144) 
(145) 

(146) 



After some rearrangement and use of trigonometric identities, we obtain from Eq. I|143|) our final result 



dP m 

dn 



2tt 



(l + 2(m|) 



0J m COS ( 



WL(w m )Y 



J m {W m ) 



tan 9 



(77l|) {qUJ m ) UJr, 
7T7 2 



cos 9 [J m {W m )f 



(147) 



in terms of the observed frequency uj m = mu>o. This distribution shows a parity violation contribution proportional 
to cos 9. The standard result for the angular distribution of the spectral power given in Ref. |45| is directly recovered 
when £ = 0. Let us observe that the expansion parameter £ appears amplified by the factor m, which turns out to be 
very large in some astrophysical settings, as we will discuss in section VIII. 



B. Total power radiated in the m i/l -harmonic 



In the previous subsection we computed the unpolarized angular distribution of the averaged power radiated into 
the m th harmonic, up to first order in £. Such linear terms vanish when we integrate on the whole solid angle, as it is 
evident by symmetry arguments. In this subsection we are interested in the leading contribution to the unpolarized 
(angular) integrated spectral averaged power distribution, which is of second order in z = £lu. We expect that such 
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second order contribution will still be amplified by factors having powers of to. In order to perform the expansion let 
us introduce the notation 



g(z) = ~r== < W 2 4 m [2m/3n(z)} - [n~ 2 (z) - (3 2 ] J dx.J 2m (x) 



in such a way that we can write 



P m = q I^ [g(z m ) + g(-z m )} = P m ,+ + Prr 

ZK 



Expanding in powers of z we have 



F lll = q -^(2g(0) + z 2 m g"(0)) 



2R 

The term g(0) is readily obtained as 

g{Q) = \^p 2 j' 2m {2mP)- 1 - 2 
The computation of g"(0) is more laborious. It can be written as 



2m/3 



dxJ2m(x) 



d 2 g(z) 




dg(n) 


! d 2 g(n) 


dz 2 


z=0 


dn 


„=1 dn2 



Thus we need to calculate the first and second derivatives of g(n) evaluated at n — 1. Our results are 

r2m/3 



dg_ 
dn 

£g 

dn 2 



= (i + F 



dx J 2m (x) + 2m/3 (l - I3 2 ) J 2m (2mf3), 



f' Zmp (2m B 

= ~2p 2 J' 2m {2m(3) - (1 + P 2 ) ^ dx J 2m (x) + i^-J- ) J' 2m {2rr, I). 



which lead to 



d 2 g(z) 



dz 2 



= 2m/3 1 ' 2 J 2m (2m/3) + {Am 2 ^ 2 - 2) p 2 J 2m (2m/3 ). 



z=0 



Going back to Eq. (|150|) giving the total power in the m th harmonic, we obtain our final result 



P m = ^j 2 /^ m (2m/?)- 7 - 



/•2m/3 

~ 2 / dxJ 2m (x) 

JO 



(148) 

(149) 
(150) 

(151) 

(152) 

(153) 
(154) 

(155) 



+ (iu m ) (3 [j- 2 [mJ 2m {2m(3) + 2m 2 0J' 2m {2m0)] - 0J' 2m {2m0)] j . (156) 

When f = we recover the well established result for the total power radiated in the m th harmonic ^5|. Again wc 
notice the dominant amplifying factor m 2 in the fourth term in the RHS of Eq. l|15(jfl . 



C. Polarization index 

Now we make an expansion to first order in £ of the polarization index in Eq. (|137J) . To this order 

W m x - (l + Afw) m/3sin6>, 

together with 

Jm(W m \) ~ J m (m/3sin0) + A£wto/3 sin# J' m (m/3 sin 9), 
J' m (W m x) J' m {mf3 sin 9) - \& 



^2 g^ n 2 q ^ 

J' m (mf3 sin 0) +m — — - — J m {mB sin 9) 



(3 sin 6 



(157) 

(158) 
(159) 
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where we have substituted the second derivative of the Bessel function arising in Eq. i|159|) via its differential equation. 
Substituting in Eq. I|137[) we obtain 



Rm — 



l-joj \ PJ' m {W m -) - cot 0J m (W m _) 
1 + |w/ PJ' m {W m+ ) + cot 9 J m (W m +) ' 

The basic building blocks can be rewritten as 

pj' m {Wm-) - cot ej rn (w m -) = P- + &Q-, 
PJ m (W m+ ) + cot J m (W m+ ) =P+- lwQ+, 

where 



P- = (3 J' m (m0 sin6) - cot dJ m (m0 sin 0), 
Q 



01 q2 ^ 

/3(1 + mcosd)J' m (m/3 sin0) + m : — J m (ra/3sin0), 



sm( 



P + = f3J' m {mfisin9) + cot 9J m (mf3 sin 9), 

Q+ 



02 q2 2_ 

(3(1 — m cos 0) J^mfl sin 6) + m : — J m (mf3 sin 6). 



1-^(2-^-% 



In terms of the above quantities, and to first order in £, we arrive at 
Here we have introduced 

R 



P + 



P- flJjn (jn/3 sin 9) — cot 9 J m (m(3 sin 9) 



P+ fiJ'm ("1/3 sin 6») + cot 9J m {m/3 sin 0) ' 



(160) 

(161) 
(162) 

(163) 

(164) 
(165) 

(166) 
(167) 
(168) 



which corresponds to the polarization index in the case £ = 0. Using Eqs. i|163!166[) we can evaluate the first order 
correction in £ to R m 



P_ 



Q + 2 cos 2 9J m (m/3 sin 0) — mf3j 2 sin 3 9J' m (mf3s\a{ 



sin^0 



[cot0J m (m/3sin0)] 2 - [(3 J r ^(m(3 sin 9)f 



(169) 



The polarization index determines the Stokes parameters, given by Eqs. (|135|l . Using Eqs. jlfiOjl . (|168|> and (|169|l . 
we get in particular 



1 + i? r 2 , 

which gives a neat signature of the Lorentz violation 



- 4R m _ ~ [/3sin0J^(m/3sin0)] 2 - [cos0J m (m/3sin0)] 2 

— 4^ — — 2£(jj j J! 

[/3 sin QJ' m (m{i sin 0)] + [cos 9J m (m/3 sin 0) J 



(170) 



VII. CONTRIBUTIONS FROM DIMENSION SIX OPERATORS 



Up to this point we have considered the Myers-Pospelov model as an exact one, which allows us to make quantitative 
predictions of the ensuing effects. The inclusion of higher dimension LIV operators in the action will certainly modify 
the results obtained for the higher order contributions in the corresponding parameters, i.e. and fj n+1 with 

n > 0. The investigation of additional dimension six operators is an interesting problem of its own, but it is beyond 
the scope of this paper and we do not consider it in any detail here. Nevertheless, it is important to emphasize that 
the contributions of these operators do not produce birefringent effects, and thus this type of effects by dimension five 
operators can be clearly disentangled from dimension six contributions by polarization measurements, for example as 
indicated in Eq. (|1T()|> . 

We can perform a rough estimation of the induced effects due to the addition of dimension six operators in the 
Maxwell action for the total power radiated into the m th harmonic, on the basis of the following considerations. In 
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the first place, since the propagation of the electromagnetic field is basically characterized by the refraction index, 
the addition of dimension six operators leads to a modification dominated by a term proportional to uj 2 



n(Xz) — > n(Xz) — n(Xz) + ^2^ 2 , 



(171) 



where £2 = &/M 2 is the parameter associated to dimension six operators, and n(Xz) is given by Eq. (|44|l . The total 
radiated power into the m th harmonic remains given by Eq. (149), but with modified g(Xz m ) functions 



g(Xz m ) — > g(Xz rn ) = g(n{Xz m )) F(Xuj m ), 



(172) 



where the real function F(Xu> m ) — F(X£,Lo m , £2^™) depends upon to only through u> m , and satisfies F(£ = £2 = 0) = 1. 
Thus we have 



P„ 



2R 



[g{n(z m ))F(uj m ) + g(n{-z m )) F(-oj m )] , 



instead of Eq. I|149ll . The corresponding expansions are now 



F(Xcj m ) = 1 + AX& m +\B? + C£ 2 ) < + . 

where we assume that A, B and C are of order one together with 

dg(n) 



g(n(z m )) = g(z m ) + £ 2 w r , 



dn 



+ 0(M- i ) 



(173) 



(174) 



(175) 



6=0 



Since we already have a factor £2 oc M 2 in front of the above derivative we can further set there £ = 0, which moans 
that it is evaluated at n = 1, and thus 



g(n{z m )) = g{z m ) 



dg{n) 
dn 



0(M~ 3 ). 



(176) 



Replacing in Eq. (|173l) and recalling that g(z) ~ g(0) + zg'(0) + ^z 2 g"(0), together with g'(0) = [dg(n)/dn] n=1 , we 
finally obtain for the total radiated power 



Pm = 



q 2 u m 

2R 



2ff(0) 



1+U +6 w 



£V(0) + 2(2^+£ 2 



dg(n) 
dn 



(177) 



The factor of g(0) in the RHS shows corrections of order ux^/M 2 that we can safely neglect. The second square 
bracket also contains contributions arising from dimension five and six operators. Using Eqs. I|151|) . (|153|) and l|155|) 
the corresponding magnitudes can be written as 



X 



m5 



2^2 



X m Q — 4:UJ 



/.2m/? 

m 2 j~ 2 J 2m (2mfi) + / dx J 2m (x) 
Jo 

2mf3 

dx J2m(x) + m~f~ 2 J 2m (2mP) 



(178) 
(179) 



Using the corresponding expressions for large to given in Eqs. l|B2f) . and the further approximation l|B6|) for the case 
m/m c < 1 we obtain 



X. 



ril6 



£2 



[l + 0.36 m 2 / 3 7" 2 ] 



X m5 e [l + 0.26 m 4 /3 7 -2] 
Thus, the effects of dimension six operators are negligible in the whole range of the spectrum if 



£2 



(180) 



/i 2 < 1. On 



the other hand, if 



6 



/£ > 7 they are significant even in the high energy limit, near the cutoff frequency where 



to = to c = 37 3 /2. At the moment there is not enough experimental evidence to set this point. For example, Ref. [3£ 



gives the boundaries 



< 10 12 and —10 6 < £2 S 10 4 - I n the following we restrict ourselves to the Myers-Pospelo\ 



model with no admixture of dimension six operators. 
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VIII. HIGH ENERGY LIMIT 

In Table [I] we present a rough estimation of the relevant parameters associated to some observed cosmological 
objects corresponding to SNR and GRB. There r [l.y] is the distance of the object to the earth, 7 is the Lorentz 
factor of the charged particles, B [Gauss] is the average magnetic field producing the synchrotron radiation, uj c [GeV] 
is the cut-off frequency and loq [GeV] is the Larmor frequency. In all cases the cut-off frequency uj c has been estimated 
from the radiation spectrum fitted by a synchrotron model in the corresponding reference. This, together with the 
magnetic field B allows us to estimate the Lorentz factor 



where M is the mass of the charged particle. From the above we further obtain the Larmor frequency 

Lo [GeV] = 0.6 x 10~ 17 f ^ ) B[Gauss]. (182) 



7M 

In the case of Crab Nebula we adopt the typical values given in Ref. [5^. For Mkn 501 we consider two possible 
models for synchrotron radiation where the emitter particles are either protons [54[ or electrons |55j | . In the latter case 
we use the radius of the orbit r' = l/u>o = 1.5 x 10 15 cm and the magnetic field to obtain the Lorentz factor. Finally 
we consider the GRB021206 (z w 1.25). According to Ref. |5(j this object has a very compact core with a radius of 
the order of 1 km and a magnetic field « 10 12 Gauss. The synchrotron emission region is about 10 8 km from the core 
[57]], so that we estimate the magnetic field to be 10 4 Gauss using the transport law Br = const.. From Ref. [5^ we 
take the cut-off frequency as oj c = 1 MeV. 

As indicated in Table [I] the radiation of interest is dominated by very high harmonics 10 15 < to < 10 30 exhibiting 
also high ratios of m/j, typically in the range 10 10 < m/7 < 10 22 . The corresponding values for 7 imply also (3 « 1. 
In this way, the high harmonics present in the synchrotron radiation of these astrophysical sources and the 7 factor of 
the corresponding radiating charges highlight the relevance of the large m and 7 limit with the constraint (TO/7) 2 3> 1 
to study Lorentz violating effects. In the following we present the most interesting results of this limit, from the 
phcnomenological point of view. 

A. Large m approximation 

We start from Eq. (|1 1 1|> by recalling the expressions for each helicity contribution to the integrated power spectrum 

P m x = q2 ^" l 1 + J (Xz \ { W 2 n 2 (\z m )J' 2m [2m(3n{Xz m )] 



[■2mn(\z m )f3 ) 

[1- (3 2 n 2 {\z m )] / dxJ 2m (x)\, (183) 







where n(Xz) are given in Eq. I|44(l . According to Ref. the important point here is the presence of two different 
regimes of what is called synergic synchrotron-Cerenkov radiation, depending upon (l — /3 2 n 2 (Az m )) ^ 0. Using 
expressions (116(1 and Q44)). this basic combination can be written in terms of the electron E and photon uj m energies 
as 

1 - (3 2 n 2 = (u 2 /E 2 ) (l + 2u4f 2 ) - 2u 2 J 2 + 2fjE ^ffE 2 - 2Au, m f [l - 2fjE - (u 2 /E 2 )] 

+ O(fi 2 + 0(f) + 0(i 3 ) + 0(f)i 2 ), (184) 

which contains the correction terms we are interested in. 
Moreover, in our case we have 

(3 2 n 2 {z) > f3 2 n 2 (-z), (185) 

so that we have three different possibilities for the total power P m = P m + + P m ~- (i) n 2 (Xz)/3 2 < 1 for A = ±1, (ii) 
n 2 (Xz)P 2 > 1 for A = ±1 and (iii) n 2 (+z)(3 2 > 1 and n 2 (— z)(3 2 < 1. They are combinations of two possible situations. 
Let us recall the corresponding expressions of the integrated power spectrum (|183|) for these basic situations. 



21 



TABLE I: Typical parameters of some relevant astrophysical objects 





r 


7 


B 




too 


m 


m/7 


CRAB 


10 4 


10 9 


10" 3 


io- 1 


10 -80 


10 29 


1Q 20 


(Mkn 501) p 


10 s 


10 11 


If) 2 


10 4 


1Q -29 


10 33 


1Q 22 


(Mkn 501) e 


10 s 


10 11 


10" 1 


10 4 


1Q -29 


10 33 


1Q 22 


GRB 021206 


10 10 


10 5 


10 4 


10 -3 


1Q -18 


10 15 


10 10 



1. Case 1 > j3 2 n 2 (Xz). 



Introducing the notation 



l-n 2 (Xz m )(i 2 



2to a < 



2/3 



together with the large to expressions given in Eq. (|B2|I of Appendix B we obtain 



q 2 uj r , 



( 3 



V3irR(l + n 2 (\z m )) \2m X c 



2/3 



m/m\ a 



dx K 5/a (x) - 2 [ — — 
5/31 ' \2fh > 



2/3 



K 2 / 3 (m/m Xc ) 



(186) 



(187) 



where the required values of n = n(\z m ) should be substituted in the sequel. Using the bremsstrahlung function 
defined in i|B4|) and the asymptotic limits given in Eqs. {B5|, (|B6() and (|B7|I . we obtain 



1/0 



1 - 2 



/ 3 



\2rh 



2/3' 



for m/fhxc >> 1, and 



P„ 



q 2 co (3) 1/6 T(2/3) 
TrR(l + n 2 ) 



2rh\ c 



Ac: 



2/3' 



-m/ihx 



(188) 



,1/3 



(189) 



for m/rh\ c « 1. These expressions are similar to the standard ones up to the factor 1/ (l + n 2 ) which might 
contribute to the corrections in the LIV parameters. Equation l|188|l exhibits rh\ c as the cut-off harmonic, which also 
depends on the LIV parameters. 



2. Case 1 < P 2 n 2 (Xz). 

Using the results of Eqs. 1B9|) in the Appendix we obtain the expression analogous to Eq. I|187fl 



P m (n) = 



q 2 mio (n 2 j3 2 — 1) 



R 



l + n 2 



A(x) - 2 (n 2 (3 2 - 1) 



Ai'(-x) 



(190) 



As shown in ^(| the expressions (|187[) and l|19U|) are continuous in n(3 = 1 . The Airy functions of negative argument 
involved in this sector are oscillating functions in the large x limit. 



B. Angular distribution of the m th harmonic to first order in £ 

To compute the unpolarized average angular distribution of the power radiated into the m th harmonic in the high 
energy limit we can start from Eq. (|147fl . By performing the required expansions we get 



dP m u 2 m q 2 



dtt 2ir 



(l -m£ui m cos^ |[j; n (m/3 sin 6»)f +cot 2 6> [J m (mf3 sin 9)] 2 j . (191) 
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This expression shows an anisotropy, which is an effect of first order in £. The emission is suppressed for m£u} m cos 9 ~ 1 
and enhanced for m^uj m cos 9 ~ — 1. Note that both effects are amplified by the presence of factor m = cj/u>q which 
can be very large in some regimes. When m is large, and since we have already expanded around £, we can use the 
standard result that the radiation becomes confined to a small angular range A9 ~ m" 1 / 3 around 9 = tt/2 0. Thus, 
if we consider the radiation in the frontiers of the beam, i.e. 9 ~ n/2 ± to -1 / 3 , the anisotropy becomes significant 
when 

Wm ~(^/^) 1/5 . (192) 



C. Total power radiated into the m th harmonic for large m to second order in £ 

The total power radiated into the m th harmonic in the large m limit can be computed directly from Eq. (|156f) . 
which gives 



p 



q 2 u„ 
R 



2(3 2 4 m (2mP) - 7 - 



2m/3 



dxJ 2m {x) + 2£ 2 



7 



^(2m/3) , 



(193) 



where we only keep the dominant term proportional to m 2 in the fourth term of the RHS in Eq. I|15tj|) ■ Substituting 
in Eq. i|193|l the asymptotic representations l|B2|) for the Bessel functions for n = 1, with m c being the corresponding 
critical value 



7 3 , 



(194) 



we get 



q 2 uj„ 



m c / m 
m \m. 



— k I — I — —K 2/3 [ — 



m 
m. 



2£ 2 



muj m (3 
7 



(195) 



For large m and (1 — 1 ) > the behavior of the power radiated in the m th harmonic can be obtained using the 
asymptotic expressions (|B5ll . (|B6|) and (|i37|) . For the case m/m c 3> 1 and 7>1 the result is 



2 / \ 1/2 

f wo / m N 



20Fi? V7 



1 - 



2£ 2 



mu) m (3 
7 



,-2^1/37^ 



(196) 



Since the regime of interest has m/7 3> 1, we keep the contribution from the LIV part that is amplified by the 
factor m 2 /7 2 , obtaining 



1 f m \ 



1/2 



2 \7T7 / 



J5 



1 - 3nE +—n 2 E 2 
4 



l + 2£ 2 



-2m/37 3 



(197) 



where we have rewritten ojq/R in terms of the electron energy E using l|15f) together with (|llf> . We also set /3 = 1 in 
all factors where no divergence arises. 

The complementary range m/m c <C 1 yields 



^9r(2/3) fq 2 B 



27. 



1+r 



~ 2 / mw m 



7 



,1/3 



(198) 



In Eqs. I|197() and (|198|l the Lorentz factor 7 is given either by Eq. or J3TJ, according to the value of fjE 3 //j, 2 



D. Polarization in the large m limit 

The general expression for the polarization index R m , from which we can calculate the reduced Stokes parameters, 
is given in Eq. H167|) . Since the radiation is mostly concentrated around the plane of the orbit, we will examine the 
different limits of 



(199) 
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In the large m approximation and in the case 1 — (3 2 > 0, one can write the Bessel function and its derivative in terms 
of the Macdonald functions as stated in Eqs. HB2|) of the Appendix. Using i|B6l) and (|_B7f) the latter can be further 
approximated according to — 1, where m c is given by Eq. (|194Jl . When m/m c 3> 1 the power of the radiation is 
exponentially damped and we have 



J m {mf3) 
J' m (m(i) 



7, (200) 



which leads to 



R m (6 = ir/2) = 1-^-. (201) 

P 7 



On the other hand, when m/m c <C 1, we have 



J m {mP) fl\ 1/3 T(l/5) ^ 1/3 



J> m (m(3) \A) r(2/3) 

yielding 



m l '\ (202) 



R m {6 - vr/2) - 1 - 4l/3/jr(2/3) t2 • (203) 

In both cases there appear amplifying factors such as m and m 4 / 3 , whose effects have to be evaluated in each separate 
situation. Let us also remark that the limit £ = reproduces a linear polarization as dictated by the value R 2 m = 1. 



E. Averaged degree of circular polarization for large m 

We can often assume that the relativistic electrons {(3 ~ 1) have an energy distribution of the type [5l| 

N{E)dE = C E- p dE, (204) 

in some energy range E\ < E < Ek, where typically p ~ 2 — 3. Let us define the circular analog of the degree of linear 
polarization introduced in Ref. |5l| as 

(p + H-p_H) 

(P+(w)+P-(a;)) 

where P-t(w) is the total power distribution per unit frequency and polarization A = ±1. The average (...) on the 
energy distribution of the charged particles is calculated as 



/>OC 

(Q(E)) =C E-p Q(E) dE. 
Jo 



(206) 



In dlllfl we have introduced the power emitted into the m th harmonic P m \ , such that the total power emitted with 
polarization A is 

Pa = V P„A = / PnA dm = / — du, (207) 

which produces 

P x (lj) = —. (208) 

We will carry the calculation only to first order in the LIV parameters. To this end we start from Eq. Ijl87(l 
and set [3=1 = n(Xz m ) everywhere, except in the critical terms involving rh\ c from where the corrections arise. 
We should also take into account that most of the radiation comes from the terms with m « rh\ c 3> 1, where 
#2/3(1) = 0.49, «(1) = 0.65. Thus the dominant term is 

PAH=P^y /3 ^f— V (209) 



\rh\c ) m \rh\c 
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where D is a constant. Starting from the definition (|205|l . making an expansion in u 
it^Xc = m\c( u )) we find 



n = & 



dPx(u>) 
du 



£=0, 77=0 



([^(^=0,^0) 



+ o(| 2 , £fj,...). 



\£ui and recalling that 



(210) 



Since n(u) = \/l + u 2 + u, we can rewrite d/du in terms of d/dn. Besides, rh\ c depends only on the combination n(3 
so that we can further go to d/d/3 obtaining 

2/3 

£ =0, fl=o \l + n 2 nj | =0i - =0 d/3 



<*Pa(w) 



du 



1 \ ' m 

K(m/m c ) 

m r J m 



(211) 



where the first parenthesis on the RHS gives 1 upon evaluation. In the second parenthesis we have already taken the 
limit n — 1, t) — 0. In this case we also have m c = 37 3 /2 and E = fi'y. Therefore, the argument of the bremsstrahlung 
function k(x) becomes 

2 LU 



-7 



A 7- 



(212) 



m 

m c 3 too 

where A 2 = 2(Mv/(3qB). The above leads to 7 = Ax" 1 / 2 . Again, a successive change of independent variables yields 

d/d/3 = (3 7 3 d/d-f = -2 A 2 d/dx, (213) 

2 /3 

where wc have set [3 — 1 in the corresponding factor. Finally, we have to analyze the term xm c ' that multiplies 



k(x) inside the square bracket of Eqs. (|209|l and 12111) . Substituting 



3772 = 3A J 



-3/2 

xmV 3 = (3A 3 /2) 2 / 3 is a constant that cancels when taking the ratio in (|210|l . Thus we are left with 

cLk(x) 



4~ 



ftuJ 



1 2 we find that 



(214) 



(k(x)) ' 

where the average (...) over E has been replaced by one over x, via the change of variables E — [i A x^ 1 ^ 2 , with all the 
constant factors cancelling in the ratio of the two integrals. It should be pointed out that the factor A 2 arising from 
the derivative in the numerator survives after taking this ratio. For the purpose of making an estimation of Tl(p) we 
take the energy range to be < E < 00. We then obtain 



Using the expression |^1| 



x M k(x) dx 



3^ 



ftuJ 



Jo 



• x (p-3)/2d?M dx 
ax 



qB J f °° X (P- 3 V 2 k{x) dx 



10 M + 2 

and comparing with Eq. I|217|) . we finally get 



r 



2 3/ V 2 3 



/!+ 1/3 > -1, 



n. ( ^ ) my.). 



with 



n(p) 



(p-3)(3p-i) (p + i)r(f + i|)r(| + A) 



p > 7/3. 



(215) 



(216) 



(217) 



(218) 



3(3p-7) (p-i)r(i + i|)r(| + ii)' 

The constraint p > 7/3 is required to avoid the divergence of the integral in the numerator of Eq. (|215|l at x = 0. The 
above result for n is the analogous for the averaged linear degree of polarization IIri, = (p + 1) / (p + 7/3) which 
is independent of the frequency j^J. On more realistic grounds, one should avoid the infinite upper limit of E, (the 
zero lower limit of x). In this case the divergence of the integral at x = 0, and hence the mathematical constraint on 
p, disappears, but the expression for II (p) becomes more complicated. 

In any case, the most important feature of the result in Eq. ((2171 is the presence of the amplifying factor (fiui/qB), 
which is independent of such details. An estimation of this factor in the zeroth-order approximation (£ = = fj), 
which is appropriate in Eq. (|217f) . yields (/Mj/qB) = w/(wo7) — ™/7j which is not necessarily a small number. 
According to the values that we have mentioned at the beginning of this section, this amplifying factor could be as 
large as m/7 ~ 10 22 . 
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IX. PHENOMENOLOGICAL CONSEQUENCES 



In our analysis of the synchrotron radiation wc have found that a careful perturbative expansion in the Lorentz 
violating parameters produced not only the naively expected factors, but also non-trivial large amplifying factors for 
some Lorentz violation effects. These unexpectedly large amplification factors open the possibility of observing such 
effects in the radiation of astrophysical sources where these factors become important. In this section we explore 
different aspects of the synchrotron radiation from astrophysical objects such as SNRs and GRBs and the constraints 
upon the Lorentz violating parameters imposed by their measurement. 



A. Radiation angular distribution 



This subsection is devoted to the study of the preferred directions together with the opening angles in the Lorentz 
violating circularly polarized synchrotron radiation. The standard Lorentz covariant results (see Ref. |45j | for example) 
are that the radiation is confined to the forward particle direction with an opening angle 69 ~ 7 — 1 . The recent bounds 
on the LIV parameters, derived from the synchrotron radiation observed from the Crab nebula, heavily rest upon 
extrapolating some of these results, based on the assumption that the LIV violations should only slightly modify the 
standard ones [25|. This has been the subject of some controversy in the literature [5!|. For this reason it is relevant 
to provide some answers arising from a specific model calculation in which the analysis of Refs. [2{| can be embedded. 

Our starting point is the angular distribution of the power emitted in the m th harmonic with polarization A given 
in Eq. I|96[l . which we rewrite here 



dP„ 



1 



4tt V /TT^ 



[Xmf3n(Xz m )J^(Wx m ) + mcotd J m {W\ m )X 



(219) 



Recalling that we are dealing with large values of m it is convenient to use the corresponding approximations in (|B2 
By doing this Eq. (|H7)|l leads to 



dP m x _ 

dn 



1 



{m\l3n(\z m ){l - (l3n(Xz m )sm6) 2 )K 2/3 - ((3n(Xz m ) sinfl) 2 ) 3 / 2 



-m cot 0(1 - (Pn{Xz m )sm9) 2 )^ 2 K 1/3 - {f3n(Xz m ) sinfl) 2 ) 3 / 2 ] } 



(220) 



Using further the approximation corresponding to m/rh c 3> 1, we can make explicit the asymptotic behavior of the 
Macdonald functions which yields 



dP„ 



dn s^ 2 v /TT4 i 

x Xf3n(Xz m )(l — (/3n(Az m ) sin 



e -2m[l-(^ri(Az m )sin0) 2 ] ' /3 

2 ) 1/4 + <»t0 {l-{pn(Xz m )sin9) 2 )- 1 / 4 



(221) 



The above corresponds to the use of the approximation (|B7|I of Appendix B. In fact this expression turns out to be 
a reasonable approximation for values of the Macdonald function argument as low as •& — 0.2 (m ~ 1). In this case it 
gives an error of the order of 15%, which is acceptable for our purposes. 

To visualize the geometry of the radiation beam we can rewrite this angular distribution in terms of the angle 
a — 7r/2 — 9. First we estimate the direction of the maximum radiated power in the m — th harmonic given by a max . 
Recalling the definition l|18tj|) for rfi\ c and introducing the variables 



2m Xc 



1/3 



sin a, 



we can write the extremum condition for (|221|) as 

3 M y ((1 + y 2 ) 2 + y (l + y 2 f 2 ) - y^\Ty~ 2 (l + V 2 ) 



1 ~ 0. 



(222) 



(223) 



We assume the above equation to be valid for /i Sj 1, for reasons to be explained a posteriori. Also, numerical 
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estimations show that y > 1 (y < 1) when \i < 1 (/i > 1), respectively. These considerations lead to the solutions 

(2to \ 2 ' 3 
~3~7 TO_1/3 ' M> 1 : (224) 

sma max — OL max — A 0.64 (2to c )~ 1/3 , /i = 1, (225) 

sina max ~ a ma2; ~ A (2m)~ 1/3 , ^ < 1. (226) 

The result (|226|l justifies the use of expression 1221|l for fj, <1. The argument d max of the corresponding Macdonald 
functions for a mnT is 



to 

u max — 



2/3 

sm 2 a maa; + ( 2 ^ ) cos 2 a,,,,,., 



3/2 

- (227) 

J m<<i 



As explained after Eq. Q221fl . we are still very close to the range where we have considered the approximation l|B7|) 
to be acceptable. 

Finally we estimate the opening angle of the radiation 56, defined with respect to the angle of maximum radiation. 
To do this we only consider the exponential factor in H221[l and determine the cutoff angle a c , defined to be the angle 
where the power decreases by a factor 1/e with respect to its maximum value. This leads to the equation 

2 / 2 \ 3/2 2 / 2 \ 3/2 

-m [1 - (/3ncosa c ) J = 1 + -m (f - (f3n cos a max ) J , (228) 

together with 

56 = 2(a c - a max ). (229) 

The radiation occurs mostly in the sector fi < 1 , thus by solving (|228H to lowest order and using the expressions l|225(l 
and iEMfl) we find 



50 = 0.98 m" 1 / 3 , 

56 c = 1.25m- 1/3 , n = l, (230) 

which effectively means 56 ~ m" 1 / 3 . 

When it is possible to ignore the photon contribution to the LIV that enters only through n(z m ), that is to say 
when we can set n = 1 as it is the case of Ref . |25j , we recover the result 

^(l-/? 2 ) 1/2 ^ 7 -\ 
for the opening angle, which was one of the starting assumptions in this reference. 

B. Birefringence 

One of the most evident features of this Lorentz violating electrodynamics is directly related with the propagation 
of the electromagnetic field and manifests itself as a vacuum birefringence effect, i.e. circular polarizations propagate 
with different velocities 



= yi+Lo 2 e±u£. (231) 

In principle this difference produces a shift in the arrival time of the different polarizations radiated by a given source, 
but in fact this time delay is very small to be observed even for a propagation in cosmological distances. A more 
sensible manifestation of this difference in the propagation velocities is the polarization of the field, which changes 
with the distance. 
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C. Radiation anisotropy 

Another interesting effect produced by the Lorentz violation is the anisotropy in the emitted radiation. According 

to Eq. Ijl92(l this anisotropy becomes apparent for oj ~ I oj 2 /£ 3 J , where the radiation is suppressed in the region 

< 9 < tt/2 {it/ 2 < 9 < tt) and enhanced in tt/2 < 9 < it (0 < 9 < it/ 2) for £ > (| < 0), respectively. This 
effect gives a clear signature of the Lorentz violation, and can be used to set a bound for £. In this respect, the most 
favorable test in laboratory conditions could be provided by an electron accelerator such as LEP. There the energy 
of the electrons is E ~ 55 GeV, the radius of the orbit is R ~ 4.25 /cm, yielding the Larmor and a cut-off frequencies 
of w - 5 x 10~ 20 GeV and w c ~5x 10~ 5 GeV respectively. The absence of asymmetry in the radiation up to the 

critical frequency should imply £ < 10 -6 GeV -1 . In fact this is a very weak bound. A more stringent one could 

be obtained from astrophysical systems such as GRBs, provided that the electromagnetic radiation of these objects 
actually corresponds to synchrotron radiation. 



D. Cutoff frequency 



The exact high energy limit for the integrated power spectrum is given in Eq. (|188|) . At this level we can identify 
the modifications induced by the Lorentz violating effects in the cutoff frequency through the term in the exponential 
e -"/ w c 5 w here 

^ c =^ (l-n 2 /3T 3/2 . (232) 

The corresponding factors expressed in terms of the electron energy, given by Eqs. (|11J) and (|184(l . are rewritten here 
to first order in the small dimensionless parameters £, r\ previously introduced 

. = M£f 1 + !,,4V (233) 



E V 2'M r 

1 - j p n *=f-2n— -2A£— . (234) 
H E 2 'M V ' 

Since the terms in (l — /? 2 n 2 ) 3 ^ 2 are already very large it is appropriate to neglect the very small correction appearing 
in ujq. In this way we have 

3\q\B E 9 ,^\- 3/2 

w c = tt^H ^7-2??— -2A£— . (235) 



2 E \E* M WAS M, 

The relation of the above result with that given in Eq. (4) of Ref. [25j can be easily found. Changing v/c by v/c(oj) 
in our definition (|18|l of the Lorentz factor we get 

7./ = {l - [v{E)/e{u)] 2 y 1/2 = (1 - n 2 (3 2 y 1/2 , (236) 
where we have introduced the index of refraction n = c/c(u). Let us rewrite 

7 . 7 = c( W ) {[cM - v(E)} [ C ( W ) + v(E)}}- 1/2 = {2 [cH - v{E)]y 1/2 (237) 
to the leading order. The dispersion relations give 

2 (cH - v{E)) = -2\^ + £ - 2^, (238) 
in such a way that Eq. (|235|l can be written as 
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which is precisely Eq. (4) of Ref. |25|. up to the polarization factor multiplying the photon frequency which has 
been subsequently taken into account in later publications. Nevertheless, in the considered energy range the photon 
contribution is negligible and the following maximum energy is obtained 



E '"'" x = ' ~5^ M 



1/3 



(240) 



It is possible to verify that 



(i-/?V) B _ 



5 -E"max 



> o, 



(241) 



so that one stays in the allowed range of radiation emission. 



E. Additional corrections 



Another factor that modifies the power spectrum and tests the dynamics of the electromagnetic field, according to 
Eq. H197fl . is 1 + 2 (j,m 2 ujo/^j • If this correction is not observed at a given frequency to wc infer that £lu 2 / (luq^/) <C 1, 
or equivalently 

i<^u Q /Lu 2 = qB/(nu; 2 ). (242) 

This bound is independent of the energy of the electrons. It is proportional to u/q/uj 2 , as in the parity violating effect 
already discussed after Eq. (|191|l . but now we have an additional factor of 7, which increases the resulting upper 
limit, thus making the boundary much less stringent. 

It is interesting to compare the bounds for £ provided by the first order anisotropic effect in (|197jl and the second 
order suppression in the power spectrum in i|191|) • According to Eq. I|192|) , the first case requires 

l< (lu 2 /lu 5 ) 1/3 = (lo/u; q ) 1/3 lo /u 2 , (243) 
while in the second case Eq. I|242ll demands 

|< 7 cj /w 2 . (244) 

Assuming that the Lorentz violation is not observed up to a frequency uj < uj c = 7 3 woi the boundary (|243|l leads to 

£ < (^M,) 1/3 lo^/lo 2 = (lj/lu c ) 1/3 ^loo/uj 2 ) < (245) 

If to ~ u> c both relations l|242|) and i|243|) give the same bound for £. 

Finally, Eq. Ijl97(l contains two additional frequency independent factors including 77— dependent corrections: 
(l + f]E 3 /2n 2 ) that comes from 7-V2 through Eq. and (l - 3f)E + 27rj 2 E 2 /4). Both test the dynamics of 

the charge. The first one dominates the linear corrections at energies that we can find in astropysical systems, and if 
they are not observed for electrons of energy E this implies 

\f]\ ^ V 2 /E 3 . (246) 



X. ACCESSIBLE PARAMETER REGIONS AND GENERAL OUTLOOK 



To close our analysis, we will determine the regions of LIV parameters that can be explored by means of the 
synchrotron radiation of the already considered astrophysical sources: Crab Nebula, Mkn 501 and GRB021206. The 
Crab Nebula is a well known synchrotron radiation emitter, while the other two are on a more conjectural level. 

Let us first consider the two LIV contributions to the power spectrum Eq. I|197|) discussed in section IX-E. The 
one containing 77 is a frequency independent factor which simply renormalizes the power spectrum. Since we do not 
have an independent determination of the density of the emitter particles we are not able to evaluate it. The second 
one, which contains £ and depends on the frequency, distorts the shape of the power spectrum thus being in principle 
measurable. If this distortion is not observed in the power spectrum we can set the bound 1242(1 . In the particular 
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TABLE II: Factors contributing to the Green function phase for the astrophysical objects 





r'/r 


LUc/Mp 


(r'/r)(co c /M P ) 


(r'/r) 2 


CRAB 


10" 6 


1Q -20 


1Q -26 


lO" 12 


Mkn 501 


HT 11 


HT 15 


1Q -26 


1Q -22 


GRB 021206 


1Q -24 


1Q -22 


1Q -46 


1Q -48 



case of the CRAB Nebulae this leads to £ < 5 x 10 -1 . Other strong sources of high energy electromagnetic radiation 
are Mkn 501 and GRB 021206. Eq. I|242|l shows that these objects could allow us to explore the region £ > 10~ 6 . 

The values of the different factors contributing to the Green funcion phase l|45l) for the astrophysical objects under 
consideration are given in Table [H] Here we take the Larmor radius as an estimation of the size of the radiating 
source r'. The last column gives the subdominant term in the phase of the Green function, (r'/r) 2 , that allows us to 
estimate if a given LIV contribution is significant for the radiation field, according to the discussion in Section III. In 
the Crab Nebulae case Table ITT1 tells us that £w = ^Lo/Mp <§; (r'/r) 2 , even for £ w 1. In this way, consistency requires 
that the phase of the Green function (|45|l reduces to u (i — n • r'). This is equivalent to set £ = 0, n = 1 in all the 
arguments of the Bessel functions. This leads to 

86 « -i-\E), uj c = qB 1 3 (E)/E, (247) 

where j(E) still includes corrections depending upon the fermionic parameter fj according to Eq. (|16[) . In others 
words, the results in Ref. p5| are recovered for this particular situation. On the other hand, the corresponding phase 
in the cases of Mkn 501 and the GRB 021206, leading to modified expressions for the cut-off frequency could, in 
principle, give information about first order contributions in £. Since these corrections are not affected by significant 
amplifying factors, their observation would require an extremely precise determination of the cut-off frequency. In 
this analysis we have considered the case of a charged particle in a circular orbit orthogonal to the magnetic field, and 
we have used our results to obtain some rough estimations of the regions where the LIV parameters can be explored. 
In realistic systems such as SNRs or GRBs, we actually have a population of charged particles with a certain energy 
distribution and different pitch angles. Thus, in order to derive more reliable bounds it is necessary to take into 
account these distributions and perform the corresponding averages. This process involves a rather detailed model 
for the astrophysical source and is beyond the scope of the present article. 

To summarize, we have presented a complete analysis of the synchrotron radiation in the rest frame of the Myers- 
Pospelov model for a charged particle in circular motion perpendicular to a constant magnetic field. The process can 
be visualized as a standard electrodynamics radiation in a dispersive parity violating media, with helicity-dependent 
refraction indices that encode the Lorentz violating corrections to the electromagnetic sector. The charged particle 
sector also presents corrections arising from the Lorentz violation, which mainly modify the dependence of the (3 
factor upon the energy of the particle. The calculation includes exact expressions for the angular distribution of 
the power H96[) and the total power pill) radiated into the m th harmonic for each polarization, together with their 
corresponding expansions to lowest order in the LIV parameter £ given in Eqs. (|147f) and i|156|) . respectively. The 
physical situations under consideration correspond to large Lorentz factors 7 together with harmonics in the range 
m ss 10 15 — 10 30 , but in such away that the ratio 771/7 is very large also. Due to this fact, a large m expansion has 
been performed in the total power radiated into the m th harmonic for the exact polarized case, Eqs. I|188I189|) . and 
in the unpolarized ^-expanded situation, Eqs. H197N1981) . The £-exact case allows us to identify the cutoff frequency 
according to Eq. i|186|) . Special emphasis has been given to the polarization of the radiation which exact Stokes 
parameters are obtained in terms of Eq. (|137f) . with their corresponding expansion written in Eq. i|167|) . We have 
also calculated the averaged degree of circular polarization for a standard electron energy distribution in the large m 
approximation in Eq. I|217|) . This quantity starts linearly in £ which, unexpectedly, is multiplied by the amplifying 
factor ([iuj)/(qB) ss 771/7. Such amplifying factors occur also in other ^-expanded quantities. Similar amplifying 
factors have been obtained in calculations of the synchrotron radiation spectra in the context of non-commutative 
electrodynamics [(J^. A study of the observational relevance of these factors is outside the scope of the present 
work, but certainly it is a matter that deserves further investigation. As in the standard case, the radiation in the 
large m limit is concentrated in the forward direction. Finally, the angles for maximum radiation together with the 
corresponding opening angles are calculated in Eqs. 12241) . <|225|) . 12261) and l|230|) . respectively. 
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APPENDIX A: CIRCULAR POLARIZATION BASIS 

Given a wave propagating in the direction k we define the ± polarization components of a vector j = {ji (to, k)} as 

1 



j- k-j k±£(kxj) , (Al) 



where we have introduced the unit vector k = k/|k|. These components satisfy the following relations 

j=J + +r, j ± -l ± = (A2) 
kxj+ = -ij+ kxj- = ij- (A3) 

with 1 being another arbitrary vector. 

Sometimes it is convenient to rewrite <|A1|1 in terms of projectors P£ defined as 



it 



p ^Jk , Pfc = - (dik - kk ± ie ijk kj) • (A4) 

The algebra of these projectors is very useful and is summarized in the relations 

P t + k = Ph > (A5) 

P&P£ = Pu> ^^ = 0, (A6) 



c 



ki il kl ' ki il 



pqr kp P qs P rt — , tpqr kp P qs P^ — ^iP st . (A7) 

Let us state that in the frame we have previously chosen we have the following basis 

k= (sin6, 0, cos ff), e ( | = (0,1,0), ej_ = (- cos0, 0,sin0), (A8) 
with the associated circular basis 

e± = —= (ey ± ie±) , (e±)* • e± = 1, (e±)* • e T = 0. (A9) 

Thus a direct calculation shows that for an arbitrary complex current j = (j X jjy,jz) the corresponding expressions 
for the currents defined in Eq. i|Al[) are 

j ± = -j= [j y ±i(j x cosd- j z sinO)] e±. (A10) 

APPENDIX B: LARGE m EXPANSION OF SOME REQUIRED BESSEL FUNCTIONS 

Using the Nicholson-Olber asymptotic representations for the Bessel functions in terms of the Macdonald 
functions K n / m , we can write, in the large m approximation [4(1 l64l|: 

1. The case 1 - n 2 /3 2 > 



)-2/3 / Q \ 2 / 3 
= (i ■ <B1) 
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1 / 3 \ 1/3 (m 
73V (2^) K2/3 (m 



J' 2m (2mnf3) ~ — ( — j K 2/3 ' "' 



/ dx.J 2m {x)~^ M1/3W, (B2) 

where the cutoff number m c signals different regimes in the behavior of the Bessel functions K v of imaginary argument. 
A useful relation among the Macdonald functions is 

2K 2 ) 3 (x) = - [K 5/3 (x) + K 1/3 (as)] . (B3) 

The bremsstrahhmg function k(z) is defined as 



k(z) = z I dxK 5/3 (x) , (B4) 



and has the following asymptotic behavior 

k(z) ~ 2 2/3 r(2/3) z 1/3 , 2 <1, 



k(z) ~ J I z 1/2 e~ z , z>l. (B5) 



We also need the asymptotic behavior of K v {z) 



K v {z)^ Y -^( Z -\\ z«l, (B6) 



2 V2 

ff v (z)~ \/|^ 1/2 e~ z , (B7) 



where the last relation is independent of za 

2. The case 1 - n 2 /3 2 < 



Here the notation is 



and we have 



x = TO 2 / 3 (n 2 /? 2 -l), (B8) 



J 2m (2m n P) - -m~ 2/3 Ai'(-x), 

J2m(t)dt = -+ / dtAi(-t) 
3 Jo 

K(x) = - 2Al '^ x) +1+1 dtAi(-t), (B9) 







x 3 

where we have introduced the analogous A(x) of the bremsstralung function k(x) 
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